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We derive tiie one-loop renormalization equations for the 
shift in the Fermi- wavevectors for one-dimensional interacting 
models with four Fermi-points (two left and two right movers) 
and two Fermi velocities v\ and 112. We find the shift to 
be proportional to (wi — V2)U'^, where U is the Hubbard-C/. 
Our results apply to the Hubbard ladder and to the ti — t2 
Hubbard model. The Fermi-sea with fewer particles tends to 
empty. The stability of a saddle point due to shifts of the 
Fermi-energy and the shift of the Fermi-wavevector at the 
Mott-Hubbard transition are discussed. 

PACS numbers: 75.30.Gw, 75.10.Jm, 78.30.-j 



Introduction - Fermi-surface properties of one- and 
two- dimensional interacting electron systems are a fas- 
cinating topic of current interest, motivated in part by 
the opening of a pseudogap in underdoped cuprates. 
Several novel effects have been found in generalized 
renormalization-group (RG) approaches to weak coupling 
models, in particulat-thc formation of quasiparticle gapsu 
and a Pomaranchul!U|and chargoJ instability in A^-patch 
2D-Hubbard models^a 

Novel Fermi-surface effects are possible also in N-ba: 
ID-dimensional interacting models, e.g. N-leg ladder 
Interest in the field of ladders systems has grown rapidly 
in the last years due to the discovery of. jSuperconduc- 
tiviy in a doped two-leg ladder materialfl'S and due to 
the fact that a rigorous weak-coupling analysis can be 
carried through. 

There have been a number of previous RG investi- 
gationS|-pf ID interacting two-band systems. Penc and 
SolyomB extended an earlier calculation of Varma and 
Zawadowskill3 and studied tiie Fermi-velocity renormal- 
ization at two-loop. Fabriziollil considered the role of the 
transverse hopping tp on the stability of the Luttinger- 
liquid (LL) state both for the case,of particles with and 
without spin. Balents and Fisherllj considered the same 
model we do, a general ID model with four Fermi points. 
There are several phases with x gapless charge and y gap- 
less spin modes (notation: CxSy) in the phase diagram 
wich were identified by Balents and Fisher considering 
the bosonized form of the leading divergent running cou- 
plings. The new feature we wish to stress is that the gen- 
eral phase diagram is influenzed by the RG-flow of the 
Fermi vectors, which we evaluate explicitly. This Fermi- 
surface effect is especially strong near a saddle point. 

It is well knownll3 that the tadpole (or Harttree) dia- 
gram contributes to the renormalization of the quadratic 
part of the original Hamiltonian, thus leading to a renor- 
malization of the chemical potential in the case of a one- 
band model. Alternatively, one might wish to work with 



a constant density of particles. This can be achieved 
by the introductiop_of appropriate counter terms to the 
chemical potentiallj. Here we would like to point out, 
that the tadpole diagram leads, in general, to different 
quadratic terms, which we denote by fii and fi2, for the 
two bands in a ID two-band model, see Fig. |l|. The 
renormalization of the chemical potential is given in the 
2-band case by suitable avarage of /xi and fj.2- The dif- 
ference of the two quadratic terms results, on the other 
hand, in a renormalization of the interchain hopping tp 
and thus alters the band structure. 

RG-equations - We will consider for a start a general 
two-band model with a linearized one-particle spectrum 
around the Fcrmi-points ki and ^2 and denote with vi 
and V2 the Fermi- velocities for the first band (Li and Ri) 
and second band {L2 and R2) respectively. Lateron we 
will speciiy our results to the case of the Hubbard ladder 
and ti — t2 Hubbard model. Since we will be working 
at one-loop, we will keep the Fermi-velocities constant 
throughout the calculation, even when the Fermi-points 
ki and k2 are allowed to shift. 

A list of possible coupling constants is given in Table 
p. The special couplings g2v, Qw, giv and gg exist only for 
the ti — t2 Hubbard model. The couplings listed come 
in two flavors, namely for scattering between parallel and 
antiparallel spins (| | and _L) for backward scattering, g-yij, 
and forward scattering, g-y2j, with 7 = 1,2, x,t and j =\\ 
, _L. Back- and forward scattering amplitudes for parallel 
spins are not independent, we choose g.yi\\ = 0. 

Spin-rotational invariance implies that the RG-flow 
leaves the combinations 
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9jcr ■— 25^1J_ 
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invariant (7 = l,2,x,t). Here the index p and a denote 
the charge and spin component respectively. 

TABLE I. Coupling constants for (a) Umklapp and (b) 
normal scattering processes (^2 > fci has been assumed). For 
the ti — t2 Hubbard model one has to interchange R2 ^ L2, 
compare Fig. (|l|) and (p[). 



label 


condition 


process'"' 


label 


process'-''' 


g2u 


k2 = Tv/2 


L2L2 — > R2R2 




L2R2 — > L2R2 


g2v 


3^2 + ki = 271 


L2L2 — > R1R2 


gi 


LiRi — > L\R\ 




3^2 — ki = 2tt 


L2L2 — > L1R2 


g^ 


L1R2 L1R2 


Qxu 


k2 + ki = n 


L2L2 — ^ RiRi 


gt 


L2R2 -> LiRi 


gtu 


k2 + ki = n 


L2L1 R1R2 


gs 


L2R1 — > LiLi 


giv 


k2 + 3ki = 27r 


L2L1 — > RiRi 




(3fci = fca) 


gin 


ki = n/2 


RiRi — > LiLi 
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The RG-equations to onarloop and momentum cut-ofF 
A are for general densitiesEHl 



2 
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9ip = /3 {g?p + Te9L) 9i<y 

92p = a (5tp + fE9D 92a = -/3.gL ~ ^ + ^agtpgta 

5a:p = — \9tp^ TE9t(j) 9xa = — 9xa ~ 29ta ~ ^9tp9ta 

9tp^ 9op9tp + ^^"f/*" = goa9tp- 



[9op 



where the couphng constants have been rescaled by 
Xj {-nivY + vi)), t = - ln(A/Ao) (Aq is the initial cut-off) 
and where we have introduced 



2vi 



2V2 



(2) 



{1/a + 1/ P — 1). The initial values of the rescaled cou- 
pling constants are 



-U 



2U 



9ip 



27r(wi -f V2) 



g^ya 



Tt{vi + V2) 
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for the Hubbard model. Inclusion of Umklapp scattering 
at special fillings leads to additional terms. They are 



9ip = • • • - a^L 
giu = -4a.gi„gip 



(4) 



for giu, wher|e-|the dots in the first equation denote the 
generic termiJ. We find 



gtp 

gs 



2ils ^ 



2ag 



ip 



gta = 
igxa /A - 
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At half filling k2 — ki = tt/2. Depending on /3, various 
coupling constants might scale to strong coupling, the 
divergence is in general of the form ~ g* /{t*~t). In Table 
H a complete list of the diverging coupling constants is 
given. The differences between the Hubbard ladder and 
the ti — 12 Hubbard model, which stem from the different 
^ometries of the two Fermi-seas, are pointed out in Table 

The Umklapp scattering g2u at ^2 = vr/2 is not active 
in the ti — t2 Hubbard model, since k2 > 7r/2 for this 
model. g2v is allowed in the ti — t2 Hubbard model and 
would leads to a CISO phase if (a) 3^2 + ki = 2^ and (b) 
13 > 2.4. It turns out, that it is not possible to satisfy 
both conditions at the same time. Inspecting the fix- 
point HaHii|ltonian and bosonizing the diverging coupling 
constantsll3, see Table one finds that the number of 
gapless charge and spin modes are C1S2 (lu,2u), CISO 
(lv,2v), ClSl (s,w) and COSO (xu,tu). 

RG of Fermi-wavevectors - The renormalization 
equations for the local chemical potentials are given by 
the tadpole diagram: 



Mi 



A(t;i + V2) (gip + gxp) 



1,2) 
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Note that the cut-off A occurs on the right-hand side of 
Eq. and that the fii are linear in the coupling con- 
stants. Let's consider now a renormalization step, we 
write 



(10) 



where the counter-term AfXct is introduced in order to 
enforce particle number conservation. 



for gs and 



g-r'p - 

giv = gi 
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25, 



tp 



gip = ■■■ - '^glv 

2gxp - 2agip 



(6) 



for giy (with 7' = x,t). The RG-equations can be solved 
numerically, we present in Fig. |2| a typical fiow diagram 
for some selected coupling constants. 

Phase diagram - Above RG-equations can be easily 
implemented either for the Hubbard ladder or for the 
ti — t2 Hubbard with respective dispersion relations 



El^''\k) 



zLtp — 2t cos(fc) 



E^^^\k) = -2iiCos(fc) - 2t2COs(2fc) 



(7) 



The reflection symmetry in the Hubbard ladder (ex- 
change of the two legs) forbids any coupling with an odd 
number of operators per band, see Tables ^ and ||. 

In Fig. ^ we present the phase-diagram for the ti — 12 
Hubbard model, for which the control-parameter /? takes 
the form 



Vi + V2 
2V2 



sin(A;i) + sin(A;2) 
2 sin(fc2) 



(8) 



TABLE II. Diverging 


coupling constants 


in all possible 


phases, k 


= max(a, /3). 






phase 


conditions 


divergent couplings 








^ (+00) 


C2S2 


K > 4.8 


(none) 


(none) 


C2S1 


4.8 > /? > 4.3 




(none) 


C2S1 


4.8 > a > 4.3 




(none) 


CISO 


4.3 > K 


gtp,gxp,gi,2cT 


gta,gi,2p, {gxcrf 


lu 




gip 


giu 


2ut 




g2p 


g2u 


lv*,2v*f 


K > 2.4 


9tp,gxp,gi,2<T 


gtcr, gi,2p, gxcr 


s* 




gxp,gip,gtp 


gs,g2p 


w* 




gxp,g2p, gtp 


gw,gip 


xu^,tu^ 




gtup 


gxu ^ gtua 






gtp,gxp,gi,2cT 


gta,gi,2p,gxa 



*Does not occur in the ladder due to transversal momen- 
tum conservation. 

^ Docs not occur in the ti — t2 model due to conditions 
on fci and k2. 

^ Subleading divergence '-^ g%al\/i* ~ 
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The counter term is compensated by subtracting the 
corresponding term from the interaction part of the 
Hamihonian. The size of A/Ltct is determined for the Hub- 
bard ladder by the condition {AkiVi — A/i°) 

Afci = -Ak2 , V2Aii°^ = -viAfi^ (11) 

(for the ti — t2 Hubbard model Aki — A/c2 etc.), which 
leads to 



= (-1)^^(AM2-Am2) 

AWct = — ^A/i2 + ^2_Aui . 

Eq. (|l|) and Eq. (|) lead then to 

k = {-iyA{g2p- 9ip) , 



(12) 



(13) 



which describes the flow of the Fermi- wavevector. Since 
the initial values for gip and g2p are identical, see Eq. (||), 
a non-trivial renormalization of the Fermi- wavevector oc- 
curs only when vi ^ V2, i.e. when a ^ (3. For the 
Luttingcr-liquid phase C2S2 the limit 



Afci = lim 

t — >oo 



dt' kit') 



(14) 



is well defined. In the other phases there is normally 
a plateau in Akiit) = dt' ki{t'), due to the decreas- 
ing cut-off A in the expression ( |l3|) for ki which can be 
considered to define an approximate Afc^, see Fig. |[ 

Let us consider the fiow of the Fermi-wavevector in 
some more detail for the Hubbard ladder, where fci < ^2 
and vi < V2. The relation a > (3 then leads generally to 
52p > <?ip and consequently to 



Afci < , Afc2 > 



(15) 



The smaller Fermi-sea tends to empty. 

In the limit i7 — > the renormalization of the coupling 
constants tends to zero in the C2S2 phase, and one can 
obtain an asymptotic rigorous expression for the total 
shift Afc = -Afci using {i = 1,2) 



9ip{t) ~ 5v(0) + gv(0)i 



9^piO) 



1 



Vl + V2 TT^V3-i 

and Eq. (|) for gip{0). With (H) and (|l|) we obtain 



Afc 



Tr'^{vi + V2) \Vl V2 
TT^ {vi + V2Y 



A 



• 



(16) 



The Fermi-surface renormalization effects become espe- 
cially strong near a saddle point, where vi <C ^2, i-e. 
a ^ /3. Eq. ( |l6| ) has been derived for the case of the 
Hubbard ladder. For the ti — ^2 Hubbard model one has 
just to interchange a <-> /3, both Fermi-points shift then 



in the same direction (to larger values), see the inset of 
Fig. I 

Discussion - In above discussion of the Fermi-surface 
renormalization we have kept the Fermi- velocities vi and 
V2- At two-loop, the Fermi- velocities are renormalized, 
an effect 0{U^). The change of 5k is, on the other hand, 
quadratic in U and leads therefore (for the ladder) to 
an additional reduction in ui, which is 0{U^). This re- 
duction in Vl I V2 would lead to a further enhancement of 
Afc in a self-consistent treatment. We do not, however, 
expand this point here further, since it would involve a 
somewhat uncontrolled mixture of high- and low-energy 
scales. 

It is of interest, however, to study the stability of the 
saddlepoint itself. This question has been studied exten- 
sively for the case of A^-natch models of the 2D Hubbard 
model. Honerkamp et al.a have found, that particles tend 
to leave the saddle-point region (also called 'hot spo 
in accordance to our result Eq. (16). Gonzalez et al. 
have found that the saddle point attracts the Fermi-level 
both from above and from below. 

We have therefore considered the case when the Fermi- 
level is exactly at the bottom of the first band in the 
Hubbard ladder. In this case the usual RG-approach has 
to be modified. We hgoie followed the approach proposed 
by Balents and FisherE3 (an expansion in the curvature) 
and find with vi = 



Av2gx 



fl2 = Av2g2p ■ 



(17) 



As the initial values of gxp and g2p are identical for the 
Hubbard model, see Eq. (||), one has to consider the 
flow of gxp and g2p, which is, due to the lack of scakL. 
invariance, explicitly cut-off dependent. One findsE£l 
g2p > 9xp, i-e. the Fermi- level is pushed below the bottom 
of the first band. The saddle point is therefore pushed 
above the Fermi-level by the RG-flow. 

Mott-Hubbard transition - A Mott-Hubbard tran- 
sition occurs in the half-filled one-dimensional ti — t2 
model as a function of interaction strength U . For 
^2 > ^i/2, that is for the case of four Fermi-points, this 
transition-occurs at finite values of U . Aebischer et al. 
estimatecO a critical Uc ~ 2.67ti for t2 = 0.7ti. The 
question we would like to ask now is: How large is the 
Fermi-surface shift Afc for U — Uc"! We cannot answer 
this question within our weak-coupling approach, but we 
can estimate the order of magnitude of Afc by assum- 
ing formula Eq. ( |l6|) to hold up to the Mott-Hubbard 
transition. We obtain 



Afc 



t2=0.7ti 
C/=2.67ti 



0.4 Ao = 0.4(7r-fc2) 



(18) 



where we have taken n — k2 as a physical relevant ini- 
tial cut-off. This estimate indicates that the Fermi- 
wavevector is shifted by a substantial fraction towards 
the Brillouin zone edge near the Mott-Hubbard transi- 
tion. We note that the real Afc might be even larger since 
we have neglected in this estimate the renormalization of 
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the Fermi-velocity which, as we have discussed further 
above, would enhance the effect. In principle it would be 
possible that the Mott-Hubbard transition coincides with 
/c2 ^ TT for U ^ Uc- Physically this would correspond to 
an activation of the one-band Umklapp scattering gi„ at 
the transition point. 

Conclusions - We have derived and discussed the RG- 
equations for the Fermi-wavevectors in one-dimensional 
models with four Fermi-points. We find that the renor- 
malization of the Fermi-wavevectors is quadratic in the 
interaction, suggesting substantial effects at strong cou- 
pling and near the Mott-Hubbard transition. We have 
pointed out, that the renormalization of the Fermi- 
wavevector is of especial importance near a saddle point, 
and that it leads to changes, Awi, in the Fermi- velocities 
quadratic in the coupling constant C/, whereas the usual 
two-loop contributions to Au^ are cubic in U . 

In addition we have discussed the phase diagram of 
the ID ti — t2 Hubbard model and found a new phase 
(ClSl,s/w). 
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Inll3 g2p — > Up at the saddle point, see Eq. (3.7) of Ref.t3. 
C. Aebischer, D. Baeriswyl and R.M. Noack, Phys. Rev. 
Lett. 86, 468 (2001). 
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FIG. 1. Dispersion relation of a Hubbard ladder with 
interchain hopping tp. The renormalization flow will lead, in 
general, to different "chemical potentials" for the two bands 

(a) , which in turn will lead to renormalized Fermi-wavevectors 

(b) . For the case of the Hubbard ladder, this process can be 
recast in a renormalization of the interladder hopping tp. 
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FIG. 2. Example of the RG-flow of some selected cou- 
pling constants for the half-filled ti — 12 Hubbard model with 
U = O.Ui, t2 = 0.7ti, a = 0.63, /3 = 2.4 (CISO phase, see Fig. 
^ corresponding to ki — 0.4187r, k2 = 0.9187r and vi — 1.66fi, 
112 ~ OASti. The initial values of the coupling constants are 
given by Eq. (^ . Also shown is the flow of the Fermi-surface 
shift Ak (in units of Aq/IOO). 
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FIG. 3. Phase diagram for the ti — t2 Hubbard model. 
The region denoted 'one band' is the ClSl LL-phase with 
one band empty, n = 0.5 dejiptes half-filhng. Note that the 
'lu' phase-hne is continuoustj. Inset: An illustration of the 
dispersion relation. 
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